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ABSTRACT 

We report results of a three dimensional, high resolution (up to 512 3 ) numerical investigation of super- 
sonic compressible magnetohydrodynamic turbulence. We consider both forced and decaying turbulence. 
The model parameters are appropriate to conditions found in Galactic molecular clouds. We find that 
the dissipation time of turbulence is of order the flow crossing time or smaller, even in the presence of 
strong magnetic fields. About half the dissipation occurs in shocks. Weak magnetic fields are amplified 
and tangled by the turbulence, while strong fields remain well ordered. 

Subject headings: MHD - turbulence - waves - ISM: kinematics and dynamics - ISM: magnetic fields 



1. INTRODUCTION 



The large linewidths of molecular species in molecu- 
lar clouds in our Galaxy imply the velocity dispersion 
a v of the gas is much larger than the sound speed c s 



(cr v 



1 - 10 km s" 



whereas 



0.2 - 0.3 km s" 1 ; 
Schlocrb (1997)| ). Al- 



e.g. Solomon ct al (1987), Heyer 

though magnetic field strengths arc difficult to measure 
the best estimates in such clouds give Alfven speeds va 
that are much larger than the sound speed but of o rder, 
or somewhat exceeding , the velocity dispersion (e.g. My- 
ers & Goodman (1988) L Crutcher et al (1993) , Goodman & 



Heilcs (1994), Crutchcr (1998)). The dynamics of this gas 



is of considerable astrophysical interest, as it may govern 
the rate and character of star formation in our galaxy. 

Two notions about the dynamics of molecular clouds 
have been particularly influential since the first CO maps 



vcrc made in the 1070a ; — (1) The turbulent motions in 
molocular cloude are thought to act as a "turbulon t pre s 



|Chan"- 



sure" to support a cloud against self-gravity (cf. 
drasekhar (1951)). This is motivated by the discrepancy 
between estimated cloud collapse times £ 3 x 10 6 yr, af- 
ter which it is presumed that most of the cloud would turn 
into stars (violating limits on the Galactic star formation 
rate), and estimated cloud lifetimes ^ 3 x 10 7 yr. (2) Su- 
personic, sub-Alfvenic turbulence is thought to persist for 
more than a cloud flow crossing time over cloud size L, 
t f (L) = L/a v = 10 7 yr x (L/10 pc)(er„/l km s -1 ) -1 , be- 
cause magnetic fields provide a cushion that reduces dis- 
sipation rates. In particular Arons & Max (1975) pro- 
posed that molecular cloud turbulence may be primarily 
in Alfvenic motions because for linear-amplitude waves no 
compressions are involved. Recently, both of these ideas 
have been called into question. This Letter describes high 
resolution, three dimensional numerical experiments de- 
signed to test the latter idea, with possible implications 
for the former. 

We evaluate the dissipation rate of supersonic, sub- 
Alfvenic turbulence in the context of an idealized numer- 
ical model. Our model is three-dimensional, compress- 



ible, ideal (no explicit resistivity, viscosity, or ambipolar 
diffusion), non-self-gravitating (future papers will discuss 
the effects of gravity), isothermal (a fair approximation 
for most of the material in molecular clouds) and has a 
uniform mass to flux ratio. It is also homogeneous and 
isotropic, insofar as it considers the evolution of turbu- 
lence in a periodic box, where there are no boundaries. 

This work builds on earlier results from our group and 
from others. Gammie & Ostriker (1996; hereafter GO) 
considered similar issues in a 1 2/3 dimensional model, 
while Ostriker, Gammie, & Stone (1998; hereafter OGS) 
considered a 2 1/2 dimensional model. Among other re- 
sults, GO found that purely Alfvenic turbulence quickly 
couples to other, compressive waves (see §4 below), and 
OGS found that independent of initial turbulence levels, 
magnetically supercritical clouds collap se gravitationally 



in 5-10 Myr, in the absence of stirring. Padoan & Nord- 



lund (1997)| studied the evolution of Mach 5 decaying 



MH D turbulence in 3D cl oud models with (3 — 2 and 0.02, 
and MacLow et al (1998) studied the evolution of Mach 5 
decaying MHD turbulence for (3 — 1 and 0.04; both groups 
concluded the dissipation time in such models is short. 
This work differs from these last two in that we consider 
driven turbulence and turbulence decaying from saturated 
initial conditions, thus avoiding transients associated with 
a particular choice of initial conditions. 

2. METHOD AND PARAMETERS 

We integrate the equations of compressible, ideal MHD 
using the ZEUS code (Stone & Norman 1992a; 1992b). 
The model is a cubic, periodic box of size L containing a 
plasma of uniform density po threaded by an initially uni- 
form magnetic field Bo = (Bo, 0, 0). The sound speed c s is 
constant in both space and time. Grid resolutions vary be- 
tween 32 3 and 512 3 . To allow a study of turbulent mixing, 
all models evolve a passive contaminant which initially fills 
a cylindrical volume in the center of the grid oriented with 
the symmetry axis parallel to Bo and with diameter and 
axial length equal to L/2. 
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We drive turbulence by adding velocity perturbations 
5v at time intervals At with AtL/c s = 0.001. Each Sv 
is an independent realization of a Gaussian random field 
with power spectrum \S~v^\ oc k 6 exp(— 8k/k p k) (we choose 
k p k = 8x(2tt/L)), subject to the constraint that V-5v = 0. 
Since the input perturbations are incompressive, this is a 
minimally dissipative way of stirring the model. The per- 
turbations are normalized so that the kinetic energy input 
rate E = const., and no net momentum is added to the 
box, J pSv = 0. 

In addition to AtL/c s and k p kL/2ir, two dimensionless 
parameters characterize our models: (3 and E. We study 
values of = 0.01 (strong field), 0.1 (moderate field), 
1.0 (weak field), and oo (pure hydrodynamics). Corre- 
sponding physical values of the magnetic field are given 

by B — 1.4/iG /T 1 / 2 (T/10 K) 1/2 (n H J10 2 cm- 3 ) 1/2 . In 

this Letter we present only models with Ej pqL 2 c z s = 10 3 , 
but comment on results drawn from other models. 

To transform dimensionless parameters to astronomi- 
cally relevant quantities, one may independently choose 
values for po, c s , and L. As an example, consider a cloud 
clump of size L = 2 pc, mean density nu 2 — 10 3 cm -3 , 
and temperature T = 10 K. Then the sound speed is 
c s w 0.2 km s . This implies the sound crossing time 
t s = L/c s ~ 10 Myr, and driving power E = 0.4L Q . For 
(3 — 0.01 the magnetic field strength is B = 44 pG and 
the Alfven speed is Va ~ 2 km s" 1 with the correspond- 
ing Alfven wave crossing time tA = L/va ~ 1 Myr. For 
the moderate (weak) field case the field strength is re- 
duced to 14//G (4.4/LtG), so that v A = 0.6(0.2) km s^ 1 and 
tA = 3(10) Myr. A typical velocity dispersion in a cloud 
with these properties is 1 km s _1 , so tf(L) = 2 Myr. 



3. RESULTS 

First consider driven turbulence models. Evolution of 
the total energy in fluctuations E = Ek + Eb = j (pv 2 /2 + 
(B 2 — B 2 )/(8n) for /3 = 1 models computed at resolutions 
of 32 3 through 512 3 shows that in each case, E rises steeply 
and then reaches a final, saturated value. The amplitude 
of Ej pL z c 2 s depends quite sensitively on the numerical res- 
olution; grids of 32 3 , 64 3 , 128 3 , 256 3 and 512 3 zones give 
saturated energy levels of 9.4, 13, 16, 17, and 18 respec- 
tively. Note there is a clear (although slow) trend towards 
convergence in these numbers, although since we have not 
used identical realizations of the forcing spectrum in each 
case we cannot measure the rate of convergence precisely. 
In the saturated state the dissipation rate balances the in- 
put power E, and one may define the dissipation timescales 
tdi SS = E/E and tf iss = Ek/E; these may be compared 
to the flow crossing time at the scale X p k = L/8 at which 
the turbulence is driven, tf(X p k) = \ p kj \J2Ek (which we 
hereafter abbreviate as tf). 1 All the models saturate at 
times ~ tf. For the (3 = 1 model at 512 3 resolution, we 
find t dlss /t f = 0.75 and t%Jt s = 0.54. 

Where does the energy go in the numerical models? One 
route to dissipation is via shocks. ZEUS uses an artificial 
viscosity to capture shocks, and so the shock dissipation 



rate can be measured by integrating the work done by ar- 
tificial viscosity over space and time. This accounts for 
about 50% of the dissipation. Another route is through a 
turbulent cascade li ke that which occurs in inco mpressible 
hydr odynamic (e.g. Landau fc Lifshitz (1987) ) and MHD 
(e.g. Goldrcich & Sridhar (1995)) turbulence; there non- 



linear interactions transfer energy to progressively smaller 
and smaller scales until a dissipation scale is reached and 
the energy is thermalized. Since the present models in- 
clude no explicit resistivity, viscosity, or ambipolar diffu- 
sion to thermalize the energy at small scales, energy is 
finally lost through numerical effects at the grid scale. 2 
A completely satisfactory study would include astrophys- 
ically appropriate values for the microscopic diffusion co- 
efficients and close the energy equation, but we have not 
done so here because the task is prohibitively expensive. 
This is one of the major challenges for future work. 

How does the dissipation rate vary with magnetic field 
strength? Figure la shows the evolution of E for various 
(3 at a resolution of 256 3 . The amplitude of E increases 
monotonically with field strength (decreasing [3)\ hence 
dissipation decreases as field strength increases. In Ta- 
ble 1, we give the values for the energy in the saturated 
state (averaged over time t = 0.2 — 0.3t s ), as well as the 
saturated-state dissipation times for the four models dis- 
played. From the values in the Table, the change in E 
with (3 is not large, amounting to only a ~ 30% increase 
in the E saturation amplitude as (3 varies from oo to 0.01. 
The dissipation times for saturated turbulence all lie in the 
range ~ 0.5 — 0.8i/, with slightly longer dissipation times 
for stronger-Bo models. 

The structure of driven compressible MHD turbulence 
changes as the field strength is varied. Figure 2 shows the 
logarithm of the density along three faces of the computa- 
tional volume, representative magnetic field lines, and an 
isosurface of the passive contaminant after saturation for 
both [3 — 0.01 and (3 = 1 models computed at a resolution 
of 256 3 . In both cases, the density is compressed into small 
scale knots and filaments; in the (3 = 0.01 model these are 
elongated in the direction parallel to the field. The mass 
(volume) weighted mean of \og(p/po) in the strong mag- 
netic field model is 0.28 (-0.29), whereas for the weak field 
model it is 0.20 (-0.22), indicating the density contrasts 
are larger for strong fields at fixed turbulent Mach num- 
ber. The maximum density in the strong field model is 83; 
for the weak field model it is 44. The passive contaminant 
is confined to a narrow range of flux tubes for (3 = 0.01, 
indicating that cross-field diffusion is small; for f3 = 1 it 
diffuses isotropically. 

There is a tendency toward equipartition of kinetic and 
magnetic energy in all the models. From Table 1, the 
turbulent magnetic energy SEb is between 30%-60% of 
Ek- In the weak field case, significant amplification of the 
magnetic field is produced by the turbulence, so that after 
saturation the energy in the fluctuations in the field is ten 
times larger than that in the mean field. In the weakly 
magnetized model the field lines are thoroughly tangled 
(Fig. 2b). In the strong field model the field lin es are 
relatively well ordered (Fig. 2a), as expected (e.g. Weiss 



tf is often referred to as the "eddy turnover time" for incompressible turbulence; for compressible flows the present terminology is preferred 
Ex pliritly resistive experiments done hy us rapture another gQSj of the dissipation The relative importance nf the different routes tr 



dissipa ;ion at small scales will depend on the precise values ot the microscopic diffusion coefficients (e.g. Biskamp fc Welter (iy9Uj ] 
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Table 1 

Dissipation Characteristics of Saturated MHD Turbulence 



model 


P 


E/pL 3 c 2 3 


E K /pL s c 2 s 


SE B 
E K 




t K 

L diss a 
*/ 


tdec a 

*/ 


t K 

dr. c. '■' 

*/ 


A 


0.01 


20.3 


13.0 


0.56 


0.83 


0.54 


0.82 


0.65 


B 


0.1 


18.9 


11.8 


0.61 


0.74 


0.46 


0.69 


0.39 


C 


1.0 


17.0 


12.9 


0.32 


0.70 


0.53 


0.58 


0.37 


D 


oo 


15.4 


15.4 





0.69 


0.69 


0.40 


0.40 



a t diss , t dis3 , tf, t dec , and t dec denned in text 



(1966)). 

Next consider models of decaying turbulence. The ini- 
tial conditions are taken from the saturated driven models 
presented above. Figure lb shows the evolution of E for 
decay from saturated initial conditions for various mag- 
netic field strengths. At late times the decay of E follows 
a power law in time, with index between .8-0.9 (consis- 
tent with the finding of MacLow et al (1998) ) . This implies 
that the dissipation time varies with time. We define decay 
times tdec (tdec) as t ne time taken for 50% of the initial en- 
ergy (kinetic energy) to be lost; values for the decay time 
in these decay runs are given in Table 1. For all models, 
the decay times are in the range 0.4-0.8 tf, comparable to 
the range of steady-state dissipation times. 

The decay rate measured here could in principle differ 
substantially from decay simulations that begin with un- 
saturated initial conditions. To investigate this possibility, 
we have computed the decay of supersonic turbulence from 
initial conditions in which the magnetic and velocity field 
perturbations are taken from the saturated, driven model 
A, but the density is reset to a uniform value. The result is 
plotted as a dashed line in Figure lb. The corresponding 
decay times are td ec /tf — 0.80 and t^ ec /tf = 0.68, nearly 
identical to those for Model A's decay. 

Finally, to make contact with other studies of decaying 
MHD turbulence, we have performed simulations which 
begin with a uniform density and magnetic field, and ve- 
locity perturbations that follow a k~ 2 spectrum normal- 
ized to have the same initial energy as our driven turbu- 
lence simulations at saturation. The result is shown as a 
dotted line in Figure lb; the decay times for this model 
are td ec /tf — 1.0 and t^ ec /tf = 0.6, again comparable to 
the other dissipation times we have found. Thus we con- 
clude that turbulent decay times are not strongly affected 
by specifics of initial conditions. The energy decay times 
found for 2.5D models (OGS) are a factor 1.5-1.75 times 
larger than those obtained here with 3D models. 



4. DISCUSSION 

Taking together the results of all the models presented 
here, our conclusion is that compressible MHD turbulence 
dissipates rapidly - in less than one flow crossing time 
at the energy-containing scale - regardless of the field 
strength, and of the details of the initial or ongoing energy 
input characteristics. For our model 2 pc cloud, the tur- 
bulent dissipation time is always less than 1/4 Myr when 



the turbulent scale is 1/8 the size of the cloud, and less 
than 2 Myr when the largest turbulence scale is the same 
size as the cloud. For GMCs, with flow crossing times 
^5—10 Myr and the largest energy-containing scale prob- 
ably comparable to the cloud scale, our results imply that 
without continual energy input, the observed nonthermal 
linewidths would decay in less than the cloud lifetimes. 
The energy input rate required to keep the turbulence go- 
ing is not large, however, amounting to only 0.4 L Q in 
mechanical power, input at a scale 0.25 pc, for the param- 
eters of our model cloud. 

Why are magnetic fields unable to reduce dissipation 
in supersonic, sub-Alfvenic turbulence? One might have 
expected that nonlinear but still incompressive Alfven 
waves could safely store a significant fraction of the en- 
ergy. But for the nonlinear amplitudes Sva ^ O.Iva and 
j3 = c 2 s /v\ < 1 conditions which are likely to obtain in 
molecular clouds, couplings between the MHD wave fam- 
ilies are strong. Even a circularly polarized Alfven wave, 
which is an exact nonlinear solution to the equations of 
motion, is d ynamically unstable to decay into compres - 
sive motions flgagdeev fc Galeev (1969)1 |Goldstcin (1978)|) . 
Thus a spectrum of Alfven waves of nonlinear amplitude 
can be quickly converted to compressive motions, which 
decay rapidly. 

The dissipation of turbulent energy may be an im- 
portant heat source within molecular clouds. From our 
driven-turbulence models, where energy is input at 1/8 
the size of the cloud, we find that the volume dissipa- 
tion rate of energy is E/L 3 = 7 .Bpoa 3 ,/ L; for clouds 
where the energy is primarily contained on the largest 
scale possible, the dissipation rate might be reduced by 
a factor 8 (for a given velocity dispersion). The corre- 
sponding volume- averaged heating rate is T tur b = 5.8 x 
lQ~ 27 (n H /l cm~ 3 ) (a v / km s -1 ) 3 (Lj pc) -1 erg cm _3 s _1 
for the smaller-scale turbulence case, with a factor 8 re- 
duction possible. Except for conditions where the velocity 
dispersion approaches the sound speed (e.g. dense cores 
of size ~ 0.1 pc, cf. Goodman et al 1998), this aver- 
age turbulent heating rate exceeds the cosmic-ray heating 
rate, and in large clouds where the velocity dispersion is 
large, Tturb ma y compete with photoelectric heating (e.g. 



Spitzer (1978)). Locally, however, turbulent dissipation 



may greatly exceed other sources of heating; future studies 
will characterize the localization of turbulent heating. 
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Fig. 1. — Total energy in fluctuations E versus time for various j3 at a numerical resolution of 256 3 Left: forced models 
Right: decaying models. Solid curves - decay from saturated state; dashed curve - decay of density-reset saturated state; 
dotted line - decay of k~ 2 velocity spectrum 
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Fig. 2. — Images of the logarithm of the density (colors) on three faces of the computational volume, representative 
magnetic field lines (dark blue lines), and isosurface of the passive contaminant (red) after saturation. Left: /3 = 0.01 
Right: (3 = 1 




- Images of the logarithm of the density (colors) on three faces of the computational volume, rep 
: field lines (dark blue lines), and ^osiirface of the passive contaminant (red) after saturation. Left. 



